function [alpha, beta] = fn_fit(v2, v4, v)

m = [v2*transpose(v2), v2*transpose(v4); ...
    v4*transpose(v2), v4*transpose(v4)];

RHS = [v2*transpose(v); v4*transpose(v)];

tv = m\RHS;
alpha = tv(1);
beta = tv(2);

end